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Recently Me Gettrick [T] introduced and studied a discrete-time 2-state quantum walk 
(QW) with a memory in one dimension. He gave an expression for the amplitude of 
the QW by path counting method. Moreover he showed that the return probability of 
the walk is more than 1/2 for any even time. In this paper, we compute the stationary 
distribution by considering the walk as a 4-state QW without memory. Our result is 
consistent with his claim. In addition, we obtain the weak limit theorem of the rescaled 
QW. This behavior is strikingly different from the corresponding classical random walk 
and the usual 2-state QW without memory as his numerical simulations suggested. 

Keywords: 4-state quantum walk, quantum walk with memory, localization 

1 Introduction 

The quantum walk (QW) is a counterpart of the classical random walk. The QW has two 
types like the random walk, that is, one is discrete-time and the other is continuous-time. 
In this paper we focus on the discrete-time case. Let Z be the set of integers. The 2-state 
QW on Z has been well studied by various authors [2l El SJ E] - For example, the weak limit 
theorem was given in [6j [7] . 

Mc Gettrick pQ introduced and investigated 2-state QWs with one-step memory (or also 
called "with 2nd order") on Z. We consider the walk as a 4-state QW without memory 
by rewriting the state \n2,ni,p) (= \ri2) ® |ni) ® \p)) with [n,2,n\,p) G I 2 x {0,1} in his 
model as |m, rii — ri2 + 1 + p) (= \ni) C3> |ni — n% + 1 + p)) with (ni, ri\ — rii + 1 + p) G Z x 
{0,1,2,3} in our 4-state model. The integers m, n% G Z correspond to the position and 
p E {0,1}, n\ — ri2 + 1 + p G {0, 1, 2, 3} mean the (coin) state or the chirality. He derived 
an expression for the amplitude of his QW by path counting method. Let P(X t — x) be the 
probability that the quantum walker, X t , exists at position x G Z at time t starting from 
the origin. He claimed that F(X2t — 0) > 1/2 for any t by induction on the time step. 
Corresponding to his result, we calculate the stationary distribution of X t and see that the 
return probability P(-X~2t = 0) converges to 2 — \[2 = 0.58578 • • • , as t — > oo, by using the 
Fourier analysis. The result is consistent with his claim. In general, if limsupj^^ P(Xj = 
0) > 0, we say that localization occurs. Therefore our result insists that localization occurs for 
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an initial state. Related with quantum physics, localization of the QW has been investigated 
8, 9, [TOl HU [T21 [13] . Moreover he numerically reported that the behaviour of the 2-state QW 
with one-step memory, equivalently the 4-state QW without memory, is remarkably different 
from the corresponding classical random walk and the ordinary 2-state QW without memory. 
By the Fourier analysis again, we show that X t /t converges weakly to a random variable as 
t — > oo. The limit measure is described by both a (^-function corresponding to localization and 
a density function. A similar limit theorem was presented in Inui et al. |10j for the 3-state 
Grover walk. 4-state models were also studied in Brun et al. [14], Venegas-Andraca et al. 
[T5] and Segawa and Konno [TB]- Their 4-state walks have a 3-direction shift operator (i.e., 
left, right, and center). On the other hand, the walker in our 4-state walk moves to the left 
or right. 

The rest of the present paper is organized as follows. Section 2 treats the definition of 
our 4-state QW. In Section 3, we present the limit theorems as our main result. Section 4 is 
devoted to proofs of the theorems. Summary is given in the last section. 

2 Definition of 4-state QW 

In this section we define a 4-state QW without memory corresponding to the 2-state QW 
with one-step memory introduced and studied by Mc Gettrick pQ. In the paper, the time 
evolution for his QW (called case (c) there) is written by the following operators: 
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where \m, ni,p) is the state, a, b, c, d € C are amplitudes and C is the set of complex numbers. 
Here the meaning of each element of |7i2,7ii,£>) is as follows; 712 € Z is the previous position 
(corresponding to the one-step memory), 711 6 Z is the current position, and p € {0,1} is 
the coin state (or chilarilty). The state space of his 2-state QW with one-step memory is 
composed of the set of the following vectors: 

\n — 1,71, 0) , |ti— 1,71, 1), |ti + l,7i,0), |n + l,7i,l) (n € Z). 

By rewriting |?i2,7ii,p) in his setting as |tii,tii — n-i + 1 + p) in our 4-state model, the state 
space of the 4-state QW is composed of 



|n,0), \n,l), \n,2), \n,3) (n G Z). 
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For the usual 2-state QW (without memory) with left chirality state |0) (= \L)) and right 
chirality state |1) (= \R)), we put, for example, 

|0)= T [1,0] |1)= T [0,1], 

where T is the transposed operator. In a similar fashion, for our 4-state QW (without mem- 
ory), we put 

|0) = T [1, 0,0,0], |1)= T [0, 1,0,0], 
|2) = T [0, 0,1,0], |3) = T [0, 0,0,1]. 

By using this notation, the operator defined by (fTJ) becomes 

a\n,0) +b\n,l) , 

c \n, 0) + d \n, 1) , ( ) 

o|»,2) + 6|n,3), 1 ' 

c \n, 2) + d \n, 3) , 

\n,2), 
K 1) , 
|n,Q), 
|n,3) , 

|n-l,0), 
|n-l,l), 
jn+1,2) , 
|n+l,3), 

where SC2 is shift operator S of his model. The amplitudes of C\ are described in Table [T] 
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Table 1. Amplitudes of Ci 
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The matrix expression of C\ becomes 

" a c " 
^ b d 

1 ~ a c ' 

b d 

In a similar way, the permutation matrix C2 determined by C2 is 

"0010" 

^ 10 

2 ~ 1 ' 

1 
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That is, we exchange |0) with \2). The matrix S corresponding to S is given by the following: 
S = l« - !) H ® (10) (01 + |1) + ("I ® (|2) (2| + |3> (3|). 

Therefore, the one-step time evolution operator is given by S{I®C2C\), where / is the infinite 
identity matrix. We should remark that 
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The operator S becomes the position shift operator. Then it determines the dynamics of our 
4-state QW. In other words, |0) and |1) correspond a left-mover and |2) and |3) correspond 
a right-mover. The states |0),|1),|2) and |3) correspond to the movement "right— > left" , 
"left— > left" , "left— bright" and "right— bright" , respectively (See Figure [lj. 
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Fig. 1. Relation between the states |0) , |1) , \2) , |3) and the walk with memory. 



From now on we define the 4-state QWs. Let \x) (x £ Z) be an infinite-component vector 
which denotes the position of a walker. Here x-th component of \x) is 1 and the others are 0. 
Let \tpt(x)) G C 4 be the amplitude of the walker at position x at time t. The 4-state QW at 
time t is expressed by 

i* t >=5>>®hMz)) (=x;k^(«)>). 
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Rewriting \n±l,n,p) in Mc Gettrick [I] as \n, =pl + l+p) in our 4-state model, his time 
evolution can be described as the following 4x4 unitary matrix: 



U(= C 2 Ci) 
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Noting that |0) and |1) correspond a left-mover and |2) and |3) correspond a right-mover, U 
is divided into P and Q as follows: 
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Therefore the evolution is determined by 

\^t+i(x)) = P \Mx + 1)> + Q \^ t {x - 1)> , 



(3) 



as in the case of usual 2-state QW without memory. This implies that P is a weight of the left 
movement and Q is the right one. The probability that the quantum walker X t is at position 
x at time t, P(X t — x), is defined by 

P(X t = x) = \\\Mx))\\ 2 , 

where || |a;) || 2 = (x\x). 

For example, we calculate P(X 2 = 0) and P(Xi = 0) under condition a = b = c = —d = 
1/V2 and |^o(0)) = T [a x , a 2 , a 3 , a 4 ], l^ofc)) = T [0, 0, 0, 0] (x ^ 0), since Mc Gettrick p] 
computed the same quantity. From ([3]), we have 

1^2(0)) =P|Vi(l)) + Q|^i(-l)) 

= PP |^ (2)> + PQ |^o(0)) + QP |Vo(0)) + QQ |^o(-2)> 
= {PQ + QP) |^o(0)) 
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If ai = a 2 = a 3 



04 = 1/2, we have 

P(X 2 = 0) 



¥(X 4 = 0) 



(4) 
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The result is consistent with his claim p], i.e., V(X2t — 0) > 1/2 for any t > 0. In addition, 
if 0:3 = 1, ol\ = ol% = ot4. = 0, we have 



IMo)) = \ 



(5) 



This is equivalent to his computation (2.13) in pQ. 

In order to obtain the limit theorems, we introduce the Fourier transform y$>t(k)) (fc G 
[— 7r, 7r)) of \ip t (x)) as follows: 

|* t (fc)) = ]T e -^|<M*)). 



By the inverse Fourier transform, we have 
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The time evolution of ^(fc)) is 



\% +1 (k)) = U(k)\* t (k)) : 



where U(k) — R(k)U and 
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From (j6|), we see that 



\Mk)) = U(kY\Mk)}- 
The probability distribution can be written as 



P(X t = x) 
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In the present paper we take the initial state as 

f T [a, p, 7, <S] (x = 0) 
l v " \ T [0, 0, 0, 0] (x?0) 

where |a| 2 + |/3| 2 + I7I 2 + \5\ 2 = 1. We should note that |# (fc)) = |^o(0)). 
3 Limit theorems for the 4-state Hadamard walk 

From now on we focus on the 4-state QW with a = b = c = — d = 1/V2. That is, 
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In this case, ([2]) becomes 
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which are considered as dynamics of the model corresponding to the 2-state (usual) Hadamard 
walk. So, this QW is called the 4-state Hadamard walk here. 

In this section we present two limit theorems for QW. Both Figures [2] and [3] depict the 
typical probability distributions of the walk. 
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Fig. 2. The probability distributions at time t = 500. 
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Fig. 3. The behavior of the probability distributions for time t by density plot. 



8 Limit theorems for quantum walks with memory 

For the 4-state Hadamard walk, we obtain the following stationary distribution for any 
initial state. 



Theorem 1: 
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where z is the complex conjugate, 3?(z) is the real part of z G C and 
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For a = /3 = 7 = S = 1/2, Mc Gettrick [I] showed P(X 2t = 0) > 1/2 for any t > by 
induction on the time step. Corresponding to his result, Theorem 1 gives rimt_j. 00 P(X2t — 
0) = 2--y/2 = 0.58578 see Figured] (a). Then localization occurs. On the other hand, if 
a = —f3 = —7 = 5 = ±1/2, then Theorem 1 also implies lim f ^ 00 P(X2t = 0) = 0, see Figure 
|U(b). Therefore, localization does not occur. This is in sharp contrast to the previous case. 
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Fig. 4. The behavior of the probability ¥(Xt = 0) for even time t. 



The weak limit measure of the rescaled (usual) 2-state Hadamard walk does not have 5- 
measure corresponding to localization. In fact, the following result was given by Konno [BJ [7] 
for any initial state \ipo(0)) — T {a,(3] : for — oo < a < b < oo, 

ton P f a < ^ < ft) =J {l - (|a| 2 - \/3\ 2 + a/3 + <3/3) z} / A -(x) etc, 

where 

7r(l - I 2 )V1 - 2x 2 V V5'^2> 

and I a (x) = 1 if x E A, Ia(x) = if x ^ A. In contrast to the above result, our weak limit 
measure has a ^-measure. 

Theorem 2: For — oo < a < b < oo, 
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(a<^- <b^j = {A<5 (a;) + (c + cue + c 2 2: 2 )/k(z)} das, 



where <5o(a0 denotes Dirac's 8- function at the origin. Here A, co,ci,C2 are determined by 
initial state \ipo(0)) = T [a,/3,7, 6] as follows: 



V2 , 1 

+V2lft(a7) - (4 - 3V2)3?(/35)} 



A =1 - — + - {(V2 - 2)(|/?| 2 + |«5| 2 ) + (2 - V2)5R((a - j)0 - 8)) 



co =~ -»(o7 + /3l), 
Cl -|<5| 2 -|/3| 2 + 5R((a- 7 )(^ + 5)), 

c 2 H/?| 2 + |7| 2 - ^ + »((a - 7 )(* " /8) + «7 + 3/35). 



As for the limit density function (co + cix + C2X 2 )fx{x), see Figure [5] In addition, we 
have A = 1/ V2 in Figure [5] (a) , while we have A = in Figure [5] (b) . 
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(a) |^o(0)) = T [l/2, 1/2, 1/2, 1/2] ( b ) |^ (0)) = T [l/2, -1/2, -1/2, 1/2] 

Fig. 5. The limit density function (co + c\x + C2X 2 )fx{x). 

As a relation between Theorems 1 and 2, the following equation can be obtained 

lim P(X t = x) = A. 

The result for a 3-state QW without memory in Inui et al. [10] is similar to Theorem 2. 
Moreover, localization of a multi-state QW was reported in [TTj . 

4 Proofs of theorems 

In this section we will prove Theorems 1 and 2 in Section 3. Our approach is based on the 
Fourier analysis given by Grimmett et al. [17] . 

4-1 Proof of Theorem 1 

The eigenvalues Xj(k) (j = 1, 2, 3, 4) of U(k) can be computed as 
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The eigenvector \vj(k)) corresponding to Xj(k) is 
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V N i( k ) Xj(k) {V2X 3 {k) + e lk } {V2Xj(k)e ik + 1} 



where Nj(k) is the normalized constant. The Fourier transform \ ^/ (k)) is expressed by \vj(k)) 
as follows: 
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By the inverse Fourier transform, 
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From a similar argument in |10) . using the Riemann-Lebesgue lemma, we see 
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Combining ©, ©, (HQ), (33) with (JT2) completes the proof. 
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4.2 Proof of Theorem 2 

First the r-th moment of X t becomes 
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where D = i(d/dk) and (t) r = t(t - 1) x • • • x (t - r + 1). Let hj(k) = DX 3 (k)/ Xj(k). Then 
we obtain 
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Therefore we get 
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7 ( _4_ i 4_ ) (a;) }> dor, (13) 



Thus the proof of Theorem 2 is completed. 
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5 Summary 

In the final section we conclude and discuss our results of the 4-state QW. Mc Gettrick Q] 
introduced and investigated a new kind of 2-state QWs with one-step memory. We showed 
that his walk becomes a 4-state QW by relabeling his notation \n2,ni,p) (e.g., \ri2,ni,p) — > 
|ni, ni — n% + 1 + £>))• Similarly, any extended version of his walk with r-step memory can 
be considered as a 2 r+1 -state QW without memory. In this paper, we obtained two limit 
theorems for the 4-state Hadamard walk corresponding to the case (c) studied in his paper. 
From Theorem 1, we found that localization occurs for an initial state. Moreover Theorem 
2 implies that X t /t converges weakly to a random variable with a S- measure as t — ¥ oo for 
any initial state. One of the interesting future problems is to obtain the limit theorems of 
the QW for a general a, b, c, d g C of U and m-state. For an m-state QW different from our 
model, Segawa and Konno |16j presented a convergence theorem for a suitable limit of t and 
m — ¥ oo . So we think that it is important to clarify the relation between them. 
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